Physical and chemical inhomogeneities inside the Vela SNR shell. 
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O I We present the results of the combined analysis of three XMM-Newton EPIC 

^ ' observations of the northern rim of the Vela SNR. The three pointings cover an 

area of > 10 pc^ (at 250 pc) behind the main shock front and we aim at studying 
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with high resolution the spatial distribution of the physical and chemical prop- 
erties of the X-ray emitting plasma on this large scale. We produce count-rate 
images and equivalent width maps of the Ne IX and Mg XI emission blends. We 
also perform a spatially resolved spectral analysis of a set of physically homo- 



O 

' geneous regions. We reveal physical and chemical inhomogeneities in the X-ray 

emitting plasma. In particular, we find large variations of the O, Ne, Mg, and Fe 
abundances. In some bright knots we also find unexpectedly enhanced Ne and 
Mg abundances, with values significantly larger than solar. Our results support 
' a possible association of a few X-ray emitting knots with previously undetected 

^ ■ residuals of stellar fragments (i. e. shrapnels) observed, in projection, inside the 

■ Vela shell. 
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C3 ; 1. Introduction 

The study of middle-aged supernova remnants (SNRs) in the X-ray band is traditionally 
considered a powerful diagnostic tool to obtain information about the dynamics of shock- 
cloud interaction and the relevant properties of the interstellar and circumstellar medium. 
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However, although the bulk of the X-ray emission in evolved remnants is associated with 
the post-shock interstellar medium (ISM), there is emerging evidence for the presence of 
X-ray emitting ejecta even in middle-aged SNRs. For example, there are claims for the pres- 
ence of chemical inhomog eneities i n the Cygnus Loop, observed in X-rays both at large 
( iMiyata fc Tsunemil Il999l . but see iLevenson et al.l |2002| ) and small (ILeahyl 120041 ) spatial 
scales. The metal abundances found in the Cygnus Loop are always lower than (or con- 
sistent with) the solar values, but the presence of such chemical inhomogeneities suggests 
that the X-ray emission cannot be completely associated with the shocked ISM. 

In the Vela SNR there are indications of X-ray emitting ejecta. lAschenbach et al. Jl995h 
identified 6 "shrapnels" (labelled Shrapnel A-F), that are X-ray emitting features (with a 
characteristic boomerang shape) protruding beyond the primary blast wave. Shrapne l A and 
Shrapnel D have been recently observed with XMM-Newton. iKatsuda fc Tsunemi J2005h 
found high O, Ne, and Mg abundances in Shrapnel D, thus confirming its association with 
a fragment of supernova ejecta. As for Shrapnel A, a significant Si overabundance has been 
observed, while O, Ne, Mg, and Fe have solar or subsolar abundances. These results seem 



to ind icate that Shrapnel A is somehow different from Shrapnel D and iKatsuda fc Tsunemi 



( 120061 ) conclude that the ejecta in Shrapnel A are strongly mixed with the swept-up ISM. 



The results obtained in the Cygnus Loop and in the Vela SNR open up the possibility 
to study the produc ts of SN nucleosynthesis in large and near remnants. Because of its 
proximity (~ 250 pc, iBocchino et al.lll999l . ICha et al.lll999l ). Vela SNR is an ideal target for 



this kind of studies, since it is possible to resolve fine structures at high spatial resolution. 
The age of this SNR has been estimated to be ~ 10^ yr, in good agreement with the 
characteristic age of the pulsar PSR B083 3-45 11200 



fll988h . 
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The northern rim of the Vela shell presents a patchy X-ray emission indicating a com- 
plex interac ti on b etweer i the blast wave sh ock and several ambient inhomogeneities. In 
Miceli et al.l (120051 ) and iMiceli et al.l (l2006bl ) (hereafter Paper I and Paper II) we studied 
a small isolated X-ray knot (Vela FilD) located near the northern shock front. The spec- 
tral analysis shows that two thermal components (T/ ~ 10® K, Tu ~ 3 x 10® K) of an 
optically-thin plasma in collisional ionization equilibrium (CIE) are present (Paper I). De- 
tailed hydrodynamic modeling (Paper II) shows that in the FilD knot the cooler component 
(the cloud "core") originates in the cloud material heated by the transmitted shock front. 
The hotter component (the cloud "corona" ) is the result of thermal conduction between the 
cloud and the hotter shocked intercloud medium. Moreover, the peculiar orientation of the 
optical filament associated with the FilD X-ray knot is naturally explained by our model as 
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a result of thermal instabilities. 

An open issue is related to the study of the metal abundances behind the shock front. 
As shown in Paper I, the abundances found in the FilD region are unusual and not expected 
according to metal depletion and grain destruction behind shock wave models, showing Ne 
slightly overabundant, Fe underabundant, and solar O. It is then interesting to study how 
these abundances are distributed at different distances from the main shock front and how 
they are related to the abundances observed in the Vela shrapnels. There is, in fact, the 
concrete possibility that additional shrapnels are observed in projection inside the shell. 

In this perspective we present the analysis of two new XMM-Newton EPIC observations 
of the northern rim of the Vela SNR. The observations (named RegNE and FilE in Fig. 
H]), together with the FilD pointing, cover an area > 10 pc^ (at 250 pc) and allow us to 
study at high resolution the spatial distribution of the physical and chemical properties of 
the plasma on this large scale. RegNE was partially visible in the FilD pointing analyzed 
in Paper I a s a previously undete cted (it was not visible in the ROSAT PSPC observation 



discussed in iBocchino et al.lll999l ) bright X-ray knot, spectrally harder than FilD. As for 
FilE, instead, the ROSAT All-Sky Survey indicates the presence of a large limb-brightened 
structure corresponding to a long (several pc) and brig ht Ha filament at a la rger projected 



distance from the main shock front than FilD (see also iBocchino et al.ll2000l ). 



The paper is organized as follows: in sect. [2] we present the new data and the data 
analysis procedure; in Sect. [3] we present the X-ray results in terms of image analysis 
(Sections I3.1H3.3I) and spatially resolved spectral analysis (Sect. 13.41) . Finally, Sect. H] 
discusses the results, while our conclusions are summarized in Sect. [51 



2. Data processing 



The XMM-Newton EPIC observations presented here consist of the FilD observation 
presented in Paper I and of two new observations: i) Observation ID 0203960101 (PI F. 
Bocchino), with pointing coordinate a (2000) = 8''36™45" and 6 (2000) = -42°11'39" (i. 
e. the RegNE region), and ii) Observation ID 0302190101 (PI M. Miceli), with pointing 
coordinate a (2000) = 8''34™15" and 6 (2000) = -42°55 W^ (i. e. the FilE region). The 



observations have b een performed with the EPIC MOS (ITurner et al.l l200ll ) cameras and 
with the EPIC pn dstriider et aPboOlh camera on 2004 October 31 (0203960101) and on 
2006 November 14 (0302190101). The location of the two pointings is indicated in Fig. [H 
while the relevant information about the two data sets are summarized in Table [H 



The data were processed using the Science Analysis System (SAS V7.0). Light curves. 
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images, and spectra, were created by selecting events with PATTERN<12 for the MOS 
cameras, PATTERN<4 for the pn camera, and FLAG=0 for both. To ehminate the con- 
tamination by soft proton flares, the data were screened by applying a count-rate limit on 
the light curves (binned at 100 s) at high energies (10 — 12 keV for MOS and 12 — 14 keV for 
the pn camera). The count-rate limit is 0.18 counts per second for the MOS cameras, and 
0.3 counts per second for the pn camera. The screened exposure times are given in Table [1] 

All the images presented here are superpositions (obtained using the EMOSAIC task) 
of the MOSl, M0S2, and pn images of the two pointings and of the FilD pointing presented 
in Paper I and are background-subtracted, vignetting-corrected, and adaptively smoothed 
(with the task AS MOOT H). The contribution of the background was derived from the high 
signal-to-noise background event files El_fmOOOO_Ml.fits, El_00fm00-M2.fits , El_O OOOfm_PN.fits 



(obtained using the medium filter) described in detail by lRead fc PonmanI (120031 ). 



The exposure and vignetting corrections were performed by dividing the count images by 
the corresponding superposed exposure maps (obtained with the task E EXP MAP). Since 
the ratio between the MOS and pn effective area Arei depends on the spectral distribution 
of the incoming photons, we have taken the thermal spectrum of the source into account 
when scaling the pn exposure maps by Arei, to make MOS-equivalent superposed count- 
rate images. In particular, since in Paper I we found that in the FilD region the emission 
below 0.5 keV is associated with the soft component (at ~ 10^ K), while above 0.5 keV the 
hotter component (at ~ 3 x 10^ K) dominates, we assume an incoming thermal spectrum 
at 10^ K for the images in the 0.3 - 0.5 keV band {Arei = 0.16) and at 3 x 10^ K for 
the images in the 0.5 - 1 keV {Arei = 0.20), 0.85 - 0.98 keV {Arei = 0.26), and 1 - 1.24 
keV {Arei = 0.31) bands. The background- and continuum-subtracted line images and the 
equivalent width maps presente d in Sect. 13.31 were produced by following the procedure 



described in iMiceli et al.l ( l2006al ) . 



Spectral analysis was performed in the energy band 0.3 — 1.5 keV using XSPEC. The 
Ancillary Response Files were generated with the SAS ARFGEN ta sk, and the event files 



were processed using the EVIGW EIGHT task flArnaud et al.ll200ll ) to correct vignetting 
effects. The background contribution was subtracted from the same region positions on the 
CCD (i. e. in the detector coordinates), in order to take into account the inhomogeneous 
response of the instruments across the field of view. Spectra were rebinned to achieve a 
signal-to-noise ratio per bin > 5cr and the fittings were performed simultaneously on both 
MOS spectra and on the pn spectrum. Since the FilE observation was performed after the 
micrometeoroid damage to the MOSl camera, we used only the pn and M0S2 spectra in 
region 6 of Sect. 13.41 because this region partially covers the damaged chip 6. We also added 
a systematic 5% error term to reflect the estimated uncertainties in the calibration of the 
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according to iLampton et al.l (Il976l ) . 



instrumental effective area (IKirschl 120071 ) . All the reported errors are at 90% confidence, 



3. The data analysis 
3.1. The region morphology 



Figure [2] shows the mosaiced count-rate image of the three observations of the northern 
rim of the Vela shell in the two bands chosen for the study of FilD in Paper I: 0.3 — 0.5 keV 
and 0.5 — 1 keV (hereafter soft band and hard band). 

The bulk of the X-ray emission of the RegNE knot is above 0.5 keV, although its 
brightest part, at the northern end of the knot, is also visible in the 0.3 — 0.5 keV band. This 
bright "head" shows a bow-shaped morphology in the 0.5 — 1 keV band and lies outside the 
field of view of the FilD observation discussed in Paper I. South of this head there is an 
elongated, sharp "tail", extending for about 15' (corresponding to ~ 1 pc, at 250 pc) towards 
the center of the remnant. North of RegNE a different faint, diffuse region is visible in the 
0.5 — 1 keV band. This region, labelled Hreg in Fig. [21 does not present significant emission 
in the soft band. 

In the FilE region we instead observe significant emission both below and above 0.5 
keV. In particular, in the Eastern side of the field of view there is a region (labelled H FilE 
in Fig. [2]), where the surface brightness in the 0.5 — 1 keV band is, on average, twice larger 
than that in the 0.3 — 0.5 keV band. The X-ray emission of this region is harder than that 
in the western part of the field of view (region S FilE in Fig. [2]), where we observe a high 
surface brightness in the soft energy band. 

Figure [3] shows a combined optical (Ha, in red) and X-ray images of RegNE and FilE 
in the soft (in green) and hard (in blue) X-ray bands. In RegNE (upper panel of Fig. [3]) 
no significant optical emission is present in the wide region corresponding to the hard X-ray 
emitting knot, at odd with the FilD knot where both soft X-ray and optical emission are 
present. A couple of optical filaments seem instead related to Hreg. The lower panel of 
Fig. E] shows the intense Ha emission of the FilE region (in red). The H FilE structure 
with hard X-ray em ission, lies clearly "be hind" the optical filament (in agreement with the 
results obtained by iBocchino et al.l l2000l in the same region). This configuration seems to 
indicate that we are observing in X-rays the plasma behind a reflected shock front. According 
to this interpretation, the reflected shock was generated by the impact of the main shock 
with a dense cloud, where the slow transmitted shock originates the optical emission. Note, 
however, that (as shown in Fig. [2]) the angular distance of this region from the border of the 
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shell is quite large (> 50', corresponding to > 3.6 pc at 250 pc) and that projection effects 
may be present. Therefore we cannot rule out the possibility that the optical filament E is 
not physically associated with the X-ray emitting knot. The stratification between optical 
and X-ray emission is not visible in the S FilE region, where we observe both optical and 
soft X-ray emission, with narrow optical filaments that lie not behind, but "inside" the X-ray 
emitting region. 



3.2. Analysis of the photon energy map 

To obtain information about the thermal structure of the plasma, we produce the MOS 
median photon energy map (MPE map, i.e. an image where each pixel holds the median 
energy of the detected MOS photons) in the 0.3 — 2 keV band. The MPE map can be 
considered as an indicator of the hardness of the spectrum and of the average temperature 
of the plasma, but it also depends on the chemical composition of the plasma and on the 
local value of Ny{ (where the interstellar absorption is higher we observe a hardening of the 
spectra). Note however, that, as shown in Sect. 13.41 A^h is quite u niform in the region of th e 



Vela shell covered by our set of observations, (in agreement with iLu fc AschenbachI |2000| ). 
The median energy map is shown in Fig. H] and has a bin size 12", so we can examine the 
different structures with high spatial resolution. The map clearly shows that very soft X-ray 
emission originates in the FilD and S FilE regions, where we observe relatively low values 
of MPE (< 500 eV). The H FilE region seems instead to be, on average, hotter than S FilE 
and FilD. We observe the maxima of the median photon energy in RegNE and Hreg, (i. 
e. the regions where the bulk of X-ray emission is in the hard band). We then expect the 
plasma to be, on average, hotter in these two regions than in all the others. We can also 
focus on the overall trend of the median photon energy on a large scale. Figure H] clearly 
shows that two well-separate "regimes" (the red dashed line marks the separation in the 
figure) are present: a hard regime at North and a softer regime from FilD to the South. 



3.3. Equivalent width images 

To study the spatial distribution of the chemical composition of the X-ray emitting 
plasma, we produce equivalent width (EW) maps. These depend hnearly on the abundances, 
but are also influenced by the temperature, column density, and ionization age. However, 
since temperature, Nh, and ionization age can be regarded as uniform (see Sect. 13. 4p . the 
EW maps can be considered as reliable indicators of the metal abundances. As shown in 
Paper I, O, Ne, and Fe lines mainly contribute to the observed X-ray emission in this region 
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of the shell, but thanks to the good statistics of the new observations we now detect also 
the presence of Mg lines. More specifically, in all spectra the K-shell line complexes of O 
VII (at 0.56 keV), O VIII (at 0.65 keV), Ne IX (at 0.92 keV), and Mg XI (at 1.35 keV) 
are clearly visible. We here focus only on the Ne IX (i. e. 0.85 — 0.98 keV energy band) 
and Mg XI (1.29 — 1.45 keV energy band) emission lines, since it is possible to evaluate the 
continuum in an adjacent band (1 — 1.24 keV) and since the equivalent width of these lines 
will be less affected by local variations of Nh than the low-energy O lines. The EW images 
were constructed by dividing the background- and continuum-subtracted line images by the 
corresponding underlying continuum. To estimate the continuum under the lines, we scale 
the continuum band adjacent to the line emission by using as a phenomenological model a 
thermal bremsstrahlung at kT = 0.2 keV (absorbed by a column density 2.3 x 10^° cm~^) for 
the continuum, plus narrow Gaussian components for the lines. We verified that this ad-hoc 
model provides a satisfying description for the spectra in all the regions presented in Sect. 
13.41 (reduced ~ 1-1 — 1-8, with ~ 150 — 200 d. o. f. ). A word of caution should be spent 
for the Ne IX EW, since this emission line emerges over a "false-continuum" generated by 
the blending of the Fe XVII L-lines (which have energies in the range ~ 0.72 — 1.1 keV). The 
Ne IX EW will not therefore depend simply on the Ne abundance, but on the ratio between 
the Ne and the Fe abundances. Note also that in the continuum energy band (1 — 1.24 
keV) the emission may be contaminated by Fe and/or Ni L lines and this may yield to an 
over-estimation of the continuum and then to an under-estimation of our EWs. 

Figure shows the Ne IX and Mg XI equivalent width maps. In both maps huge 
inhomogeneities in the equivalent width are visible. At odd with the median photon energy 
map, we do not observe a North- South anisotropy, while an East- West anisotropy is present. 
In particular, three knots (corresponding to Hreg, RegNE, and partially H FilE) with 
enhanced values of EW are present in both maps. Since the EW depends linearly on the 
abundance, the presence of inhomogeneities in the equivalent width maps suggests that the 
X-ray emission in this part of the Vela shell cannot be associated with a chemically uniform 
medium. 

Figure [6] includes the spectra extracted from two regions with high and low Ne IX and 
Mg XI equivalent width. The areas of the two regions have been chosen so as to have almost 
the same counts in the continuum band. Figure E] clearly shows that the Ne IX and Mg XI 
emission lines are locally highly enhanced. 



- 8 - 



3.4. Spatially resolved spectral analysis 

To check whether the spatial variations in the EWs are due to inhomogeneities in the 
chemical composition of the X-ray emitting plasma or to other effects (e. g. variations in 
temperature and/or in N^), we perform a spatially resolved spectral analysis. As explained 
in Paper I, the poor statistics of the FilD observation did not allow us to perform an 
accurate study of the metal abundances and we aim at exploiting the higher statistics of the 
new observations to address this issue. We analyze the spectra extracted from the 9 regions 
shown in Fig. [5l The shape and size of the spectral regions were chosen in order to extract 
the spectra from chemically and physically uniform regions. To this end, we select regions 
with fairly uniform values of the Ne IX and Mg XI equivalent width (as shown in Fig. 
and, at the same time, very small fluctuations of the median photon energy (< 3%). 

Region 1 corresponds to the low-surface-brightness region immediately behind the bor- 
der of the shell, region 2 to the hard X-ray emitting Hreg knot, and region 3 corresponds 
to the brightest part (the "head") of RegNE. Regions 4 and 5 are located in the H FilE 
structure, regions 8 and 9 in the bright, soft X-ray emitting S FilE knot, while regions 6 
and 7 are located between the two parts of FilE, where we observe large values of the Mg 
XI equivalent width. We observe large EW both for Ne IX and Mg XI in regions 2, 3, and 
4. In region 5 we have large Ne IX EW and low Mg XI EW, while we observe relatively low 
values both for neon and magnesium in regions 1, 8, and 9. 

All the extracted spectra are described well by two MEKAL coraponents of an optically- 



thin thermal plasm a in collisional ionization equilibrium ( iMewe et al.lll985l . iMewe et al.lll986 



Liedahl et al.lll995l ) and the temperatures of the two components are quite similar to those 
found in Paper I for the FilD knot (i. e. T/ ~ 10^ K and Tu ~ 3 x 10^ K, respectively). 
Figure [7] shows a representative spectrum (extracted from region 4 of Fig. [5]) with its best- 
fit model and residual. As shown in Fig. [TJ strong O VIII, Ne IX, Mg XI, and Fe XVII 
line (which have energies in the range ~ 0.72 — 1.1 keV) complexes are associated with 
the hotter component. Therefore, we leave the O, Ne, Mg, and Fe abundances free in the 
hotter component. Since the line emission mainly originates in the hotter component, it is 
not possible to determine the chemical abundances in the cold component. Therefore th e 



abundances of the cooler component are frozen to the solar values (lAnders &: Grevesselll989l ). 
This assumption is in agreement with the results obtained in Paper II, where we showed that 
in the FilD knot the cold component is associated with interstellar cloud material heated 
by the transmitted shock front. Except for region 3, all spectra are described significantly 
better by this model than by a single temperature model (either in CIE or in non-equilibrium 
of ionization). 

In region 3 (i. e. the head of RegNE) one component is enough. In Paper I we showed 
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that also the spectra extracted from the tail of RegNE (regions 14 and 15 in Paper I) are well 
described by a single thermal component in CIE. However, we found (Paper I, Sect. 4.6) that 
also a non-equilibrium of ionization (NEI) model at higher temperature (~ 7 x 10^ K) can be 
appropriate to describe the spectra. We try to overcome this ambiguity by modeling the pn 
and MOS spectra of the head of RegNE (where the statistics are significantly better than in 
the tail) with a single MEKAL model and with the XSPEC NEI model. We do not obtain 
significant improvement in the quality of the fits by using the NEI model {Xcie ~ 264.8, 
with 243 d. o. f. and x%ei = 264.2, with 242 d. o. f.). We then adopt the CIE spectral 
model. 

We find an entanglement between the best-fit values of A^^^- and the temperature of 
the cooler component, Tj, with large error bars for both the parameters. To overcome this 
problem we determine the Nh value in region 3 (where only the hotter component is present), 
thus finding 7V^^^^ = 2.3+?-^ x 10^° cm~2. We then fix A^^ to this value for the fittings of the 
spectra extracted from the other regions. We checked that the assumption of a uniform 
in the whole observed part of the shell is valid. In fact, the minima obtained by letting the 
Nh parameter free are not significantly lower (according to the F-test) than those obtained 
by fixing it to the best-fit value obtained in region 

The best-fit parameters for all 9 spectral regions are given in Table [2j The temperatures 
of the two X-ray emitting components are quite similar to those found in the FilD region 
of Paper I, with nearly uniform values both for Tj ~ 0.9 — 1.2 x 10^ K and Tjj ~ 2 — 3 x 10^ 
K. We instead observe large variations (more than one order of magnitude) in the emission 
measures. This result confirms that, as shown in Paper I and Paper II, neither component 
can originate in the intercloud medium. The emission measure per unit area (EM) of the 
cold component shows low values (< 10^^ cm^^) at North, in the RegNE pointing, and 
large values (many 10^*^ cm"^) in the FilE region. The ratio EMj/EMu ranges between 
EMj/EMjj ~ 2.5, in the regions with the hardest X-ray emission (i. e. region 2, or Hreg in 
Fig. [21 and regions 4 and 5, or HFilE in Fig. [2]) and EMj / EMu > 5 in the softest regions 
(regions 7 and 8, i. e. the brightest part of SFilE in Fig. [2]). These results allow us to 
understand the double-regime in the MPE map: the region North of FilD presents larger 
values of the median photon energy because in that region the cold component is fractionally 
small. 



^For example, in region 1, the best-fit value of Njj is consistent with N^^^ and we obtain very similar 
vales of by letting the Nh parameter free {xjree ~ 222.4, with 206 d. o. f.) and by fixing it to N"^^^^ 
iXfrozen = 222.5, with 207 d. o. f.). Even in region 8, where we obtain a very low best-fit value of Nh 
(consistent with and lower than 1 x 10^^ cm^^, with Xfree ~ 259.9 and 212 d. o. f.), the quality of the fit 
does not change significantly by putting Nh — N^^^ ixjrozen — 263.0, with 213 d. o. f.). 
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The results of the spatially resolved spectral analysis (Table [2]) further confirm that the 
X-ray emitting medium is not chemically homogeneous. We observe large inhomogeneities 
in the O, Ne, Mg, and Fe abundances, with significant overabundances both for Ne and 
Mg. We remark that the metal abundances do not change significantly if we adopt the NEI 
model. As explained in Sect. 13. 3[ we expect the Ne IX equivalent width map to depend on 
the Ne/Fe abundances. We indeed find a very good agreement between the Ne IX equivalent 
width map shown Fig. [5] and the best-fit Ne/Fe abundances shown in Table [2l For example, 
we find the highest best-fit values of Ne/Fe (Ne/Fe> 5) in regions 3, 4, and 5, where the Ne 
IX EW map has three maxima, while we find Ne/Fe < 2 in regions 1 and 9, where the Ne IX 
equivalent width is lower. The agreement between the Mg abundance and the Mg XI EW 
map of Fig. [5] is very good. We find significantly enhanced Mg abundances in regions 2, 3, 
4, 6, while the Mg abundances is consistent with the solar value where the Mg XI equivalent 
width is lower (regions 1, 5, 8, 9, see Fig. Oand Table [2]). 



4. Discussion 

The combined analysis of the three XMM-Newton observations of the northern rim of 
the Vela shell reveals the presence of physical and chemical inhomogeneities in the X-ray 
emitting medium. 

Although all spectra are described well by two thermal components at ~ 10^ K and 
~ 3 X 10^ K, a wide region at < 6 x 10^^ cm from the border of the shell presents a higher 
median photon energy than that observed at larger distances from the main shock front. This 
is due to the lower values of the emission measure of the cold component and indicates that in 
this region the plasma is, on average, hotter and, if we assume pressure equilibrium, also more 
rarefied. At larger distances we instead observe bright optical filaments and large, dense. 
X-ray emitting clouds (with larger values of the emission measure of the cold component), 
together with a soft diffuse and faint emission maybe associated with small, unresolved, and 
dense cloudlets. 

Both the EW maps and the spatially resolved spectral analysis unequivocally show that 
the chemical composition of the X-ray emitting plasma is not homogeneous. These results 
clearly rule out the possibility that the whole X-ray emission originates in the shocked inter- 
stellar medium. Moreover the spatially resolved spectral analysis confirms the detection of 
several knots with overabundant Ne and Mg. This strongly suggests we are indeed observing 
ejecta. It would also be possible that we were observing clumps of wind residuals from the 
progenitor star shocked by the main shock front. However, stel lar winds from massive s tars 
typically show overabundances of light elements, like He and N (lEsteban fc Vilchezlll992l and 
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Chu et al.ll2003l ) and it would be difficult to explain overabundances of Ne and Mg according 
to this scenario. We then suggest we are indeed observing SN ejecta. This indication is 
further supported by the fact that the regions with high EWs in Fig. [S] (and high Ne and 
Mg abundances) are visually aligned toward the center of the shell. 

The detection of X-ray emitting e jecta associated with the evolved Vela SNR is not 
totally unexpected. lAschenbach et al.l (119951 ) have observed 6 shrapnels of X-ray emitting 
material outside th e Vela shell (Shrapnels A- F). Two of these have been investigated i n 
detail (Shrapnel D, Katsuda fc Tsunemi 20051 and Shrapnel A, Katsuda fc Tsunemi 2006 ). 
In Shrapnel D large overabundance s of O (0/0(^ ~ 6.4 ) , Ne (Ne/Ner;^ ~ 13.7), and Mg 
(Mg/Mg0 ~ 14) have been found by Katsuda fc Tsunemil (2005), who a l so rev ealed a slight 
Fe overabundance (Fe/Fe© ~ 1.4). As shown in iKatsuda fc Tsunemil (120061 ). Shrapnel A 
is instead very rich in Si (Si/Si© ~ 2.6), while O, Ne, Mg, and Fe have solar or subsolar 
abundances. Note that the angular distance from the center of the remnant to Shrapnel A 
is bigger by ~ 20% than that to Shrapnel D. The abundances in Shrapnel A indicate an 
efficient mixing between the fast-moving shrapnel and the swept-up interstellar medium. In 
both cases, however, the chemical composition suggests that the shrapnels are fast ejecta. 

We remark that in our case we would be observing the ffist shrapnels ever detected 
within the Vela shell (Shrapnels A-F are all outside the border of the remnant). Note, 
however, that although the projected position of the ejecta is behind the main shock front, 
it is clearly possible for these "new" shrapnels to be physically outside the border of the 
shell. Only shrapnels ejected approximately on the plane of the sky have been detected so 
far, but it is very likely that many other shrapnels have been ejected also in other directions. 
We then expect to find fragments of ejecta in other regions of the remnant. When shrapnels 
are observed in projection within the shell, the X-ray emission is not fully associated with 
the ejecta, because (since the plasma is optically thin) there is also the contribution of 
the shocked ISM intercepted along the line of sight. Therefore, the enhanced Ne and Mg 
abundances we find in regions 2, 3, 4, 6, and 7 are indeed weighted averages of the ejecta 
and ISM abundances. Nevertheless, we can compa re our abundances with those found in 
Shrapnel D (Fig. [1]) by iKatsuda fc Tsunemil (120051 ). In particular, we consider the relative 
abundances between heavy elements, which are less dependent on the spectral model than the 
absolute abundances. In regions 2 [Hreg) and especially in region 3 (RegNE) the relative 
abundances (normalized to the Fe abundance) are quite similar to those observed in Shrapnel 
D: 0:Ne:Mg:Fe ~ 2.0:5.0:6.3:1 (region 2) and 0:Ne:Mg:F e ^ 2.0:6.7:9.0:1 (region 3) , to be 
compared with 0:Ne:Mg:Fe ~ 4.6:9.8:10.0:1 (Shrapnel D. IKatsuda fc Tsunemil l2005l ). This 
result supports the association of Hreg and RegNE with "new" shrapnels. Although in 
region 4 the relative abundances do not decrease significantly (0:Ne:Mg:Fe ~ 1.5:5.4:5.8:1), 
at South (specially in regions 6 and 7), where we expect to intercept more shocked ISM 
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along the line of sight, the relative abundances become lower (0:Ne:Mg:Fe ~ 0.7:2.2:3.5:1 in 
region 6; 0:Ne:Mg:Fe ~ 0.8:2.1:4.7:1 in region 7). Note that the largest deviations between 
our shrapnels and Shrapnel D occur in the 0/Fe abundances. The cause may be that since 
the X-ray emission in the Vela SNR is quite soft, we expect the contribution of the post- 
shock ISM to be stronger in the low-energy line complexes (i. e. VII and O VII I ) than 



in the Ne IX and Mg XI lines. As for the Si-rich Shrapnel A, iKatsuda fc Tsunemil (120061 ) 
found subsolar or solar values for O, Ne, Mg, and Fe, with lower values of the abundance 
ratios (0:Ne:Mg:Fe ~ 0.4:1.0:0.8:1). Although we cannot estimate the Si abundances in 
our observations, we can conclude that our new shrapnels are intrinsically different from 
Shrapnel A. 

In the regions where the emission seems to be completely associated with the interstellar 
medium (regions 8 and 9) we find indications for a metal depletion in the X-ray emitting 
plasma, since the O, Ne, and Fe abundances are significantly lower than solar (see Table [2]). 



5. Summary and conclusions 

The combined analysis of the three XMM-Newton observations of the northern part of 
the Vela shell allowed us to obtain a high-resolution description of a large area (> 10 pc^) 
of the remnant. 

We observe hard X-ray emitting knots at North, while there are large, soft X-ray emitting 
knots associated with bright optical filaments closer to the center of the shell. 

All spectra are described well with two thermal components whose temperatures are 
almost uniform in the whole region. The observed inhomogeneities of the emission measures 
confirm that, as in the case of FilD, the two components cannot be associated with the 
uniform intercloud medium. 

We find inhomogeneous patterns of abundances of O, Ne, Mg, and Fe and we detect 
different knots {RegNE, Hreg and part of FilE) with significant overabundances of Ne 
and Mg. The observed pattern of abundances can hardly be reconnected with the only 
interstellar medium and/or clumps of wind residuals. We conclude that these knots are 
previously undetected shrapnels of ejecta which appear, in projection, to be inside the Vela 
shell. We expect the X-ray emission from these new shrapnels to be "contaminated" by the 
X-ray emission from the shocked ISM intercepted along the line of sight. Nevertheless, in the 
regions where we expect this effect to be less considerable (Re gNE and Hreg), the relati ve 



abundances are quite similar to those found in Shrapnel D by iKatsuda &: Tsunemil (120051 ) 
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In conclusion, we confirm the presence of two pliases in tlie X-ray emitting plasma (at 
~ 10^ K and < 3 x 10^ K), as found in Paper I. In addition, we argue that the metal 
abundances in the hot component are not uniformly distributed. Regions with oversolar 
Ne and Mg are reasonably associated with shrapnels of ejecta similar to Shrapnel D. On 
the other hand, the other regions are probably associated with the evaporation of shocked 
interstellar clouds, according to the scenario described in Paper II. 
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Table 1. Relevant information about the data. 



OBSJD 


CAMERA 




Mode 


Filter 


0203960101 


MOSl 


52.1/26.3 


Full Frame 


medium 


0203960101 


M0S2 


52.1/27.3 


Pull Frame 


medium 


0203960101 


pn 


48.6/18.2 


Ext. Full Frame 


medium 


0302190101 


MOSl 


39.8/32.7 


Full Frame 


medium 


0302190101 


M0S2 


40.0/28.5 


Full Frame 


medium 


0302190101 


pn 


35.8/23.7 


Ext. Full Frame 


medium 



* Unscreened/ Screened exposure time. 



Table 2. Best-fit parameters for tlie spectra extracted in tlie 9 regions sliown in Fig. [5l described witli two tliermal 
components in collisional ionization equilibrium (only one component in region 3). The parameter is fixed to 
2.3 X 10^'' cm~^ (as explained in Sect. 13. 4p . All errors are at the 90% confidence level. 



Parameter Region 1 Region 2 Region 3 Region 4 Region 5 Region 6 Region 7 Region 8 Region 9 



Ti (10^ K) 


1 06+°-°^ 

-'-■'-'"-0.02 


1.12 ±0.07 




-"-•"-"^-o.os 


1 07+003 

-'-•'-" -0.05 




97+0-04 

'^■^ ' -0.02 


1.06 ±0.02 


1 03+0-05 

-'-•"-''-'-0.04 


EM*i 


1 4+0.2 
-■-•^-0.4 


0.77 ±0.08 




4 1+1 

^•-■--0.3 


q q+0.5 

'^•-^-0.3 




6.5tr 


8.3 ±0.2 


c; c:+0.3 
O.O_0.4 


Tn (10^ K) 


9 1 +0.3 
^•-■--0.2 


2.941H' 


2.58 ±0.08 


2.621H' 


2.64+°:^ 


9 1 9+0.07 
^•-■-^-0.06 


2 29+0-1 

^•^^-0.06 




9 qc;+0.2 
^•'-'^-0.09 


EM*j 


20+°-^^ 


0.33 ±0.02 


-'-•'-'-0.2 


-'-•^'-'-0.05 


1 «4+0.15 

-"-•"^-o.o? 


o 4+0.7 
'^•^-0.2 


1 3+0-1 

-'-•'-'-0.2 


1.481^2^ 


2 4+0-2 

^•^-0.4 


0/0,, 


' ^-0.08 


1.1 ±0.1 


1.0 ±0.1 


'-'• '^-0.05 


57+0-04 
' -0.08 


40+0-01 

'-'•^"-'-0.04 


S5+0-02 

"-'•'-'^-0.03 


0.25 ±0.02 


0.25 ±0.02 


Ne/NcQ 




2.8 ±0.3 


q q+0.4 
'-'•'-'-0.8 


o 7+0.2 
' -0.1 


1 7+0.3 
-■-• ' -0.1 


1.3 ±0.1 


9+0-1 
u.y_o 2 


0.5 ±0.1 


0.5 ±0.1 


Mg/MgQ 


2±2 


q c:+0.5 

"-••^-o.g 


4.4 ± 1.3 


2.9 ±0.8 


1 q+0.6 
-'-•'-'-0.5 


9 1 +0.9 
^•-1-0.7 


2 0+0-4 


0.7 ±0.4 


5+°-^ 


Fe/Fee 


1 4+0-5 
-■-•^-0.2 


56+°-^ 




50+0 05 


n q9+0.08 
"-'•'-'•^-0.05 


0-6^°:? 


4S+0-07 

'-'•^'J-0.08 


21+0-03 

"-'•'^-'--0.02 


25+0-03 


xVd. o. f. 


222.5/207 


309.9/306 


264.8/243 


296.2/291 


337.2/263 


174.9/193 


309.0/234 


263.0/213 


281.4/218 



*Emission measure per unit area in units of lO^*^ cm 
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9:00:00.0 55:00.0 50:00.0 45:00.0 40:00.0 35:00.0 8:30:00.0 



Fig. 1.— ROSAT All-Sky-Survey image of the Vela SNR in the 0.1 - 2.4 keV band. The 
color scale ranges between and 33 counts, the bin-size is 45" and the image has been 
smoothed by a gaussian distribution with a = 67.5". The red circles indicate the XMM- 
Newton EPIC fields of view of the three observations (pointing the regions named RegNE, 
FilD, and FilE) discussed in this paper. The position of Shrapnel D is also indicated. 
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Fig. 2. — Left panel: mosaiced count-rate images (MOS-equivalent counts per second per 
bin) in the 0.3 — 0.5 keV band (bin size= 8") of the northern rim of the vela SNR. The image 
is adaptively smoothed to a signal-to-noise ratio 30. The 4 x 10^*^ cm scale has been obtained 
assuming a distance D = 250 pc. Right panel: same as left panel in the 0.5 — 1 keV energy 
band. The red and green boxes bound the field of view of the upper and lower panels of Fig. 
ini respectively. 



Fig. 3. — Upper panel: combined optical and X-ray image of RegNE. The Ha emission 
(from AAO/UKST H-alpha survey) is indicated in red. We have superimposed, in green, 
the X-ray emission in the 0.3 — 0.5 keV band (Fig. [21 left panel) and, in blue, the X-ray 
emission in the 0.5 — 1 keV band (Fig. [2l right panel). The color bars have linear scales 
and range between 25% and 60% of the maximum count-rate for the X-ray maps, while the 
optical image is scaled according to the IRAF zscale algorithm. The field of view of this 
image is indicated by the red box in Fig. O Lower panel: same as upper panel for FilE. 
The field of view of this image is indicated by the green box in Fig. [2l 
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Fig. 4. — MOS median photon energy map in the 0.3 — 2 keV band (bin size= 12"). In 
each pixel the local median photon energy of the photons detected by the two MOS cameras 
is reported. The pixels where the median is computed with less than 5 photons have been 
masked out. The color bar has a linear scale and ranges between 415 eV and 750 eV. The 
red dashed line indicates the border between the two regimes discussed in Sect. 13.21 Two 
point-like, bright X-ray sources have been masked out with black circles in the upper two 
fields of view. 




Fig. 5. — Upper panel: equivalent width maps of the Ne IX hne emission in the 0.85 — 0.98 
keV band. The bin-size is 20" and the equivalent width ranges between and 300 eV. The 
points where the continuum-subtracted line emission is consistent with the background have 
been masked out. The 9 regions selected for spatially resolved spectral analysis are also 
shown. Lower panel: same as upper panel for the Mg XI line emission in the 1.29 — 1.45 
keV band, the bin-size is 1' and the equivalent width ranges between and 210 eV. 
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Fig. 6. — pn spectrum (in red) extracted from regions with high Ne IX and Mg XI equivalent 
width (red regions shown in the EW maps in the insets) compared with the pn spectrum 
(in black) extracted from regions with low EW (black regions). This shows the differences 
between metal rich regions and ISM regions. 



OVII 




0.5 1 
channel energy (keV) 



Fig. 7. — pn spectrum (upper, in green) and M0S1,2 (lower, in black and red) spectra of 
region 4 of Fig. [5] with the corresponding best-fit model and residuals obtained with two 
thermal components. The contribution of the two components and the ions that mainly 
contribute to the line emission are indicated. 



